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ABSTRACT 

DWIYANTI, F. G., KAMIYA, K. & HARADA, K. 2014. Phylogeographic structure of the commercially important 
tropical tree species Dryobalanops aromatica Gaertn. F. (Dipterocarpaceae) revealed by microsatellite markers. 
Reinwardtia 14(1): 4 3-51. — Dryobalanops aromatica Gaertn. F. (Kapur) is an economically important timber 
species in Southeast Asia that can serve as a good model for studying the impact of the Pleistocene glaciations on the 
genetic diversity and distribution of species in tropical regions. Seven polymorphic microsatellite markers were 
analyzed in five natural populations of D. aromatica (TV = 120 individuals): Gunung Panti in Malay Peninsula, Lingga 
Island in Lingga Archipelago, Lambir Hills National Park, Limbang and Similajau National Park in Borneo. The level 
of gene diversity (H E ) for the five populations was relatively high with a range from 0.571 (Similajau) to 0.729 
( Gunung Panti). The high genetic diversity in the present study could be attributed to the larger refugia population sizes 
of D. aromatica than that of other species. The population genetic structure revealed two distinct groups: the Malay 
Peninsula-Lingga Archipelago and Borneo. This pattern suggests that populations in each geographical area might be 
the consequence of post-glacial expansion from one or two refugia, but that gene flow between different glacial refugia 
was fairly restricted. 

keywords: Dryobalanops aromatica, microsatellite, phylogeography, genetic diversity, genetic structure. 

ABSTRAK 

DWIYANTI, F. G., KAMIYA, K. & HARADA, K. 2014. Struktur filogeografi pohon komersial penting Dryobalanops 
aromatica Gaertn. F. (Dipterocarpaceae) dengan marka mikrosatelit. Reinwardtia 14(1): 43 - 51. — Dryobalanops aro¬ 
matica Gaertn. F. (Kapur) merupakan jenis pohon penghasil kayu yang bemilai ekonomi penting di Asia Tenggara, 
>ang mana dapat dijadikan model yang baik untuk mempelajari dampak Pleitosen glasiasi pada keragaman genetik dan 
distribusi spesies di daerah tropik. Tujuh penanda mikrosatelit dianalisis pada lima populasi alami D. aromatica fN = 
120 individu): Gunung Panti di Semenanjung Malaysia, Pulau Lingga di Kepulauan Lingga, Taman Nasional Bukit 
Lambir, Limbang and Taman Nasional Similajau di Borneo. Tingkat keragaman genetik (HE) untuk kelima populasi 
relatif tinggi, dengan selang dari 0.571 (Similajau) hingga 0.729 (Gunung Panti). Keragaman genetik yang tinggi pada 
studi ini dapat dihubungkan dengan ukuran populasi refugia yang lebih besar pada D. aromatica daripada spesies lain. 
Struktur genetika populasi menunjukkan dua kelompok yang berbeda: Semenanjung Malaysia-Kepulauan Lingga dan 
Borneo. Pola ini menunjukkan bahwa populasi pada masing-masing area geografik mungkin merupakan konsekuensi 
dari ekspansi postglasial dari satu atau dua refugia, namun aliran gen antara refugia yang satu dengan yang lainnya cu- 
kup terbatas. 

Kata kunci: Dryobalanops aromatica, mikrosatelit, filogeografi, keragaman genetik, struktur genetik. 

INTRODUCTION spot of biodiversity (Sodhi et al., 2004). It is 

renowned for its rich variety of plant species espe¬ 
cially in lowland rainforests (Myers et al, 2000; 
Mittermeier et al, 2004; Sodhi et al, 2004 and 
Ohtani et al, 2013). The Pleistocene glacial 


The geographic area containing the Malay 
Peninsula, Sumatra, Java and Borneo is known as 
Sundaland and represents a globally important hot 
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periods have been regarded as major factors influ¬ 
encing the geographical distribution, demographic 
dynamics, and pattern of genetic diversity of spe¬ 
cies (Comes & Kadereit, 1998; Haffer & Prance, 
2001; Hewitt, 2004; Soltis et al., 2006; Stewart et 
al, 2010; Ramirez-Barahona & Eguiarte, 2013). 
However, little is known about the impact of 
Pleistocene glaciations on the genetic diversity 
and distribution of species in tropical and sub¬ 
tropical regions, especially in Sundaland area, 
even though the impact of glaciations in temperate 
zones has been the main focus of many phylo- 
geographic studies (Cruzan & Templeton, 2000; 
Hewitt, 2000; Emerson & Hewitt, 2005). 

Dryobalanops aromatic a Gaertn. F., locally 
known as Kapur, can function as an ideal model 
species for studying the impact of glaciations in 
tropical regions, as the distribution of this species 
covers almost all of Sundaland. The species 
occurs naturally in the lowland mixed dipterocarp 
forests of Malaya, Sumatra, the Riau Archipelago 
and Borneo (Wood & Meijer, 1964; Ashton, 1964; 
Ashton, 1968; Symington, 2004), and is particu¬ 
larly abundant on deep humid yellow sandy soils 
with a propensity for ridges, at an altitudes below 
400 m (Ashton, 2004). In the Malay Peninsula, 
the species is confined to the eastern side south of 
5°N latitude e.g., Johor, Pahang, Selangor and 
Terengganu (Symington, 2004; Chua et al, 2010). 
In Sumatra, the species has been recorded from 
North West Sumatra (Angkola, Sibolga, Kelasan 
and Upper Singkil), East Sumatra (Bengkalis and 
Siak), Mursala, Lingga and Singkep Island 
(Ashton, 1982), while in Borneo it occurs in 
Sabah (Beaufort, Beluran, Kinabatangan, Labuan, 
Labuk Sugut, Papar, Penampang and Sipitang 
districts), Sarawak (Bintulu, Kapit, Lawas, 
Limbang, Marudi, Miri, Mukah and Tatau 
districts) and Brunei (Ashton, 2004; Chua et al, 
2010). The timber and the camphor harvested 
from this tree also give it economic value. The 
timber of the species is a moderately heavy and 
durable construction timber, while the camphor 
obtained from the wood attracted early Arab 
traders, at one time being worth more than its 
weight in gold in the Middle East and is used for 
incense and perfume (Ashton, 2004). However, 
the species has since 1998 been classified as 
'critically endangered' (CR Alcd+2cd, Bl+2c D 
c2.3) by the IUCN Red List of Threatened Species 
due to logging and camphor exploitation. Studies 
associated with this species are therefore of 
interest on several levels. 

One of the methods available for inferring 
genetic diversification and biogeographical history 


over the past several million years is analysis of 
molecular data of extant species based on nuclear 
microsatellites. Microsatellites as a class featurea 
variable number of tandem repeats with a repeat 
unit of new nucleotides and have been used as 
genetic markers. Because of their high mutation 
rates, microsatellite markers are especially useful 
in detecting recent evolutionary events such as_ 
expansion and migration after the last glacial 
maximum (LGM) about 20,000 yearsago 
(Jeffreys et al, 1985). 

Phylogeographical studies using microsatellite 
on tropical tree species from Sundaland area are 
very scarce and limited, instead many studies 
using nuclear genes and chloroplast intergenic 
spacers or introns, such as in Shorea curtisii 
(Kamiya et al, 2012), Shorea parvifolia (Iwanaga 
et al, 2012), and Shorea leprosula (Ohtani et al- 
2013). Those studies detected two major geneti¬ 
cally different groups of Sundaland populations: 
Sumatra-Malay and Borneo and suggested the 
presence of a scarcely forested land bridge on the 
Sunda Shelf during glacial periods in the Pleisto¬ 
cene and predominance of tropical lowland rain¬ 
forest at least in Sumatra and eastern Borneo. 
Based on these results, we wanted to test whether 
this pattern also occurs in D. aromatica popula¬ 
tions. Therefore, the objective of this study was to 
estimate the level of genetic diversity and to reveal 
phylogeographical relationships in five D. aro¬ 
matica populations using microsatellite markers. 

MATERIALS AND METHODS 
Plant materials 

Leaf tissue of adult trees or saplings was 
collected from five natural populations represent¬ 
ing three geographical regions covering most of 
the natural distribution range of the species except 
for the Sumatra mainland: Gunung Panti (1.8344° 
N, 103.9005° E) in Malay Peninsula, Lingga Is¬ 
land (0.1607° S, 104.5855° E) in Lingga Archipel¬ 
ago, Lambir Hills National Park (4.1983° N, 
114.0427° E), Limbang (4.7544° N, 115.0075° E) 
and Similajau National Park (3.5172° N, 
113.3414° E) in Borneo (Fig. 1). A total of 120 
individuals of D. aromatica was used in this study. 
Plant materials were dried in the field with silica 
gel and stored in a freezer at -80°C. They were 
subsequently used for DNA extraction. 

Microsatellite genotyping 

Silica gel-dried leaves were ground to a fine 
powder using a TissueLyser II (QIAGEN Japan, 
Tokyo). Total genomic DNA was extracted from 
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Table 1. Primer information for seven microsatellite loci used in this study 


Locus 

Primer sequence (5’-3’) 

Size (bp) 

Annealing temperature (°C) 

Dral87 

F: TCTCTCTTATCCAACTCTCTCA 

R: AGGGAACTAAAGCAGACATCAC 

102-160 

57 

m 

Dra266 

F: AGACTTAATAATGGAGGACGAG 

R: CCACAATTAGCCACCATCTTAC 

115-165 

54 


Dra426 

F: CCAACGCTGCTCAAAGTTCGTG 

R: GCTGGCTGGCATAATATAATCC 

224-244 

57 


Dra428 

F: CTATTGATGCCCTTATAGCTTT 

R: ACGAGCCTCTCTACTCTATAAT 

265-303 

57 


Dra471 

F: TCTCAGTCTCACAATCTATCCA 

R: TTTCTGTGTCATTTTAGCAACC 

220-248 

54 


Dra519 

F: TCAAGCCAGAAGAGATAGAGAC 

R: AT ATTCCTTT CAT ATT ATT GGG 

185-219 

57 


Dra569 

F: GTAAAACCAATACACGTACATA 

R: ATGGAAGTCATTTCATCTATTT 

221-247 

54 



each sampled tree using the modified CTAB 
method (Doyle & Doyle, 1990). Seven microsatel¬ 
lite loci (Dral87, Dra428, Dra426, Dra519, 
Dra266, Dra471 and Dra569) developed for D. 
aromatica by Nanami et al. (2007) were used in 
this study. The forward primer of each marker was 
labeled with 6-FAM, VIC, NED, or PET phos- 
phoramidite (Applied Biosystems Japan, Tokyo). 
Details of the markers used in this study are shown 
in Table 1. 

A Type-it Microsatellite PCR kit (QIAGEN) 
was used to amplify microsatellite loci. Multiplex 
PCR amplification was performed in a volume of 


5 pi, containing 1 x Type-it Multiplex PCR Master 
mix, 0.2 (pM of forward and reverse primers and 
about 40 ng genomic DNA. An Applied Biosys¬ 
tems 2720 thermal cycler (Applied Biosystems) 
was used under the following conditions: initial 
denaturation at 95°C for 5 min, then 31 cycles of 
denaturation at 95°C for 30 s, annealing for 1 min 
30 s and extension at 72°C for 30 s, followed by a 
final incubation at 60°C for 30 min. Annealing 
temperatures were 57°C for Dral87, Dra428, 
Dra426 and Dra519, and 54°C for Dra266, Dra471 
and Dra569. Fragment sizes were determined 
using an ABI PRISM™ 310 Genetic Analyzer 
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Table 2. Genetic diversity of microsatellite loci in five populations of D. aromatica 


Island name 

Population 

N 

K 

N e 

Ho 

H e 

F 

Malay Peninsula 

Gunung Panti 

20 

6.86 

4.27 

0.521 

0.729 

0.289*** 

Lingga Archipelago 

Lingga Island 

13 

7.14 

3.72 

0.681 

0.689 

0.007ns 

Borneo 

Similajau 

32 

5.29 

2.62 

0.616 

0.571 

-0.084*** 


Limbang 

27 

8.43 

3.52 

0.503 

0.575 

0.192*** 


Lambir 

28 

7.86 

3.38 

0.587 

0.623 

0.032*** 

Total 


120 






Grand Mean 


24 

7.11 

3.50 

0.582 

0.637 

0.087 


Number of samples (N), mean number of alleles per population (N a ), mean number of effective alleles (N e ), mean 
observed heterozygosity (Ho), Nei's mean expected heterozygosity (HE), fixation index (F) with the significance 
level for deviation from the Hardy-Weinberg Equilibrium (HWE): P < 0.001 (***),.not significant (ns). 


Table 3. Pairwise F ST values for five populations of D. aromatica with the significance level of population differentia¬ 
tion after Bonferroni correction 


Population 

Lambir 

Limbang 

Similajau 

Gunung Panti 

Lingga Island 

Lambir 

Limbang 

Similajau 

Gunung Panti 

-0.01289 ns 

0.09347*** 

0.22536 *** 

0.12556*** 

0.25598 *** 

0.21946 *** 



Lingga Island 

0.17874*** 

0.20128*** 

0.21273*** 

0.01111ns 

- 


P < 0.001 (***), not significant (ns). 


(Applied Biosystems) and visualized using Gene- 
Mapper 3.0 software (Applied Biosystems). 

Data analysis 

Basic statistics of genetic diversity, including 
number of alleles per locus (N a ), observed 
heterozygosity (H 0 ), expected heterozygosity (H E ) 
and fixation index (F), were calculated using 
GenAlEx software version 6.41 (Peakall & 
Smouse, 2006). Deviations from Hardy-Weinberg 
(HW) equilibrium were tested for each locus in 
each population using GENEPOP 4.0.10 
(Raymond & Rousset, 1995; Rousset, 2008). Pair¬ 
wise F S T was calculated to determine the level of 
population differentiation using the program AR- 
LEQUIN Version 3.5.1.2 (Excoffier & Lischer, 
2010 ). 

An analysis of molecular variance (AMOVA) 
was carried out using ARLEQUIN version 3.5.1.2 
(Excoffier & Lischer, 2010) to apportion variance 


within populations, among populations within 
groups, and among groups. Populations were 
grouped into three island groups for the analysis, 
i.e. Malay Peninsula, Lingga Archipelago and 
Borneo (Table 2). Genetic distances between 
populations were estimated using the Nei genetic 
distance, D (Nei et ah, 1983), as implemented in 
GenAlex version 6.41 (Peakall & Smouse, 2006). 
The resulting distance matrix was used to 
construct a neighbor-joining (NJ) phenogram us¬ 
ing MEGA Version 3.1 (Kumar et al, 2004). 

STRUCTURE version 2.3.4 (Pritchard et al, 
2000) was used to estimate the number of geneti¬ 
cally homogeneous populations (K) using a Bayes¬ 
ian model-based clustering method. A burn-in of 
50,000 iterations was performed followed by 
500,000 iterations of MCMC (Markov chain 
Monte Carlo) simulation. The model was run for a 
range of K-values from 1 to 7, with ten replica¬ 
tions each. The best K-value supported by the data 
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mas assessed according to the recommendations of 
Evanno et al. (2005), whereby the statistic AK was 
calculated based on the rate of change in the log 
likelihood of the data between successive K-values 
using STRUCTURE HARVESTER (Earl & von 
Holdt,2012). 

RESULTS 

Genetic diversity and fixation index 

The mean number of alleles per locus (NJ 
across seven microsatellite loci ranged from 5.29 
in Similajau to 8.43 in Limbang (Table 2). 
Observed heterozygosity (HJ was highest in 
Lingga Island with a value of 0.681 and lowest in 
Limbang with a value of 0.503. Expected 
heterozygosity (H E ) ranged from 0.571 in 
Similajau to 0.729 in Gunung Panti with a mean 
value of 0.637(Table2). 

The fixation index [F) ranged from F = -0.084 
in the Similajau population to F = 0.289 in the 
Gunung Panti population. Significant departures 
from Hardy-Weinberg equilibrium were observed 
in all the populations, except the Lingga Island 
population. F was significantly negative in the 
Similajau population, but it was positive in the 
Gunung Panti population. The negative fixation 
index value for the Similajau population was 
significantly smaller than zero and showed an 
excess of heterozygotes, while the high fixation 
index value in Gunung Panti was significantly 
larger than zero and showed an excess of homozy¬ 
gotes (Table 2). 

Genetic structure 

Pairwise F S T values (Table 3) ranged from 
-0.01289 (Limbang vs Lambir) to 0.21946 
(Gunung Panti vs Similajau). Significant genetic 
differentiation was found in 8 out of 10 population 
pairs after a Bonferroni correction was applied. No 
genetic differentiations were observed between 
Lingga Island and Gunung Panti and between 
Limbang and Lambir. 

The AMOVA revealed that 62.97% of the total 
molecular variance was within individuals (Table 
4). 4.44% was attributed to among-populations 
differences within groups and 16.84% to among- 
groups differences. The remainder (15.75%) 
applied among individuals within populations. A 
NJ phenogram (Fig. 2) shows that relationships 
among the five populations agree well with 
geographic distribution. 

The Structure analysis revealed the presence of 
distinct clustering in the five D. aromatica popula¬ 
tions. The value of AK was highest when K = 3. 


The populations can therefore be characterized by 
the following three clusters: the Limbang and 
Lambir populations mainly composed of cluster 1, 
the Similajau population mainly composed of 
cluster 2, and the Gunung Panti and Lingga Island 
populations mainly made up cluster 3 (Fig. 3). 

DISCUSSION 

Population genetic diversity 

The levels of genetic diversity (H E = 0.571- 
0.729) estimated for D. aromatica in this study 
were higher than for other dipterocarp species in 
Asian rain forests such as Shorea javanica 
(Rachmat et al., 2012) and S. leprosula (Ohtani et 
al., 2013), but comparable with S. lumutensis (Lee 
et al., 2006), Vateriopsis seychellarum (Finger et 
al, 2012), and S. obtusa (Senakun et al., 2011). At 
face value, the genetic diversity value in Bornean 
populations was slightly lower than that of the 
Malay Peninsula and Lingga Archipelago. Simi¬ 
larly, the diversity value based on microsatellite 
markers (H E = 0.709) reported by Lim et al. 
(2002) in D. aromatica in the Malay Peninsula 
region was higher than that reported by Harata et 
al. (2012) in Lambir, Borneo (H E = 0.589). The 
high genetic diversity in the present study could be 
attributed to the larger refugia population sizes of 
D. aromatica than that of other species. High 
levels of genetic diversity may also be associated 
with population history, reproductive strategy and 
life history of the species, which have maintained 
a larger effective population size (Lee et ah, 
2000 ; Lim et al, 2002). 

Three of five populations showed a significantly 
positive inbreeding coefficients (Table 2). The 
positive values could be attributed in part to the 
presence of null alleles, or inbreeding between 
relatives (biparental inbreeding), restricted seed 
and moderate pollen dispersal, or a Wahlund 
effect due to two or more breeding subunits inside 
a given population (Papi et al, 2012). Lim et al. 
(2002) reported that the high fixation indices 
detected in D. aromatica populations in the Malay 
Peninsula, could be explained by a Wahlund 
effect. It is possible that all populations in the 
Peninsular Malaysia originated from one large 
population, which later sub-divided into several 
sub-populations. Since no evidence of Wahlund 
effect and no calculation of the null allele were in 
the present study, biparental inbreeding in 
D. aromatica could be regarded as a possible 
reason for the high fixation index in Gunung Panti 
population (Malay Peninsula). 
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Fig. 2. Neighbor-joining (NJ) phenogram based on pairwise F S t among the investigated D. aromatica populations 



Borneo 


(A) 


Malay Lingga 
Peninsula Archipelago 


(B) 


Fig. 3. Bayesian clustering analysis of five populations ofD. aromatica using STRUCTURE (K= 3). (A) Bar plot for 
each individual showing three clusters: cluster 1 composed of Lambir and Limbang (black), cluster 2 composed 
of Similajau (gray) and cluster 3 composed of Gunung Panti and Lingga Island (white). (B) Tree Plot showing a 
representation of the genetic distance among those 3 clusters. " 1" represents the black bars (cluster 1), "2" gray 
bars (cluster 2) and " 3" white bars (cluster 3). 


Phylogeographic structure 

The phylogeographic structure of D. aromatica 
was estimated by calculating the level of popula¬ 
tion differentiation. It was found that the level of 
genetic differentiation increased with geographic 
distances. Significant pairwise values were mostly 
found for populations from different islands. For 
example, a low level of genetic differentiation was 
detected between the Malay Peninsula and Lingga 
Archipelago populations and also among popula¬ 
tions within the Borneo region, while a high level 
of genetic differentiation was detected between the 
Malay Peninsula/Lingga Archipelago population 
and Borneo populations. At the Last Glacial maxi¬ 


mum (LGM; approximately 0.02 MYA), Penin¬ 
sular Malaysia, Sumatra, Java and Borneo Islands 
were connected by the exposed Sunda Shelf 
(Voris, 2000) and at this period, some studies 
suggested a savanna corridor across the Sunda 
shelf separated rainforest of Sumatra-Malay area 
and Borneo (Kaars, 1991; Kaars & Dam, 1995; 
Bird et al., 2005). Thus, the present study sug¬ 
gests that populations in each geographical area 
might be the consequence of post-glacial expan¬ 
sion from one or two refugia, but that gene flow 
between different glacial refugia was fairly re¬ 
stricted. This suggestion is also supported by 
Iwanaga et al. (2012) who reported a very low 
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migration rate in S. parvifolia population after 
divergence of the Sumatra-Malay and Borneo 
imply that contacts of the Sumatra-Malay and 
Borneo population were limited. Therefore, when 
Sunda Shelf was exposed, the expanding savanna 
must have been strong but imperfect barrier to 
gene flow between Sumatra-Malay and Borneo. 

The genetic differentiation patterns observed in 
the present study are also evident in the result of 
the STRUCTURE analysis and the neighbour¬ 
joining (NJ) phenogram. Both results show 
different patterns of D. aromatica genetic structure 
at the cluster level. STRUCTURE analysis divides 
the structure of D. aromatica populations into 
three clusters: Limbang-Lambir, Similajau and 
Gunung Panti-Lingga Island, while the NJ tree 
shows two distinct groups: the Malay Peninsula- 
Lingga Archipelago and Borneo populations. At 
the subcluster level however, both approaches 
show the same patterns: Lambir-Limbang, Simila¬ 
jau and Gunung Panti-Lingga Island (Fig. 2). 

The clear distinction between the Malay 
Peninsula and Borneo populations observed in the 
present study was also found in S. curtisii by 
Kamiya et al. (2012), where it was reported that 
long-term population persistence with limited seed 
dispersal and low levels of gene flow may act as a 
driving force of divergence between the Malay 
Peninsula and Borneo populations. A Study of 
population structure in S. parvifolia (Iwanaga et 
al, 2012) also suggests that the clear separation of 
the two groups was due to limited genomic 
exchanges between the two groups and that 
reproductive isolation has developed between 
them. In a study of phylogeography of 
S. leprosula, Ohtani et al. (2013) reported the 
presence of a longitudinal genetic cline both 
within the Malay/ Sumatra population group and 
the Borneo group, indicating that the genetic 
structures in Sumatra and the Malay Peninsula are 
not independent from each other due to long-term 
land connections between Sumatra and the Malay 
Peninsula across the shallow Strait of Malacca 
during glacial periods (current average depth about 
25 m). Our result gives an additional support for 
the inference of non-forested land bridge at LGM 
between Sumatra-Malay Peninsula and Borneo. 

CONCLUSION 

The high genetic diversity in the present study 
could be attributed to the larger refugia population 
sizes of D. aromatica than that of other species. 
The phylogeographic structure revealed two 
distinct groups: the Malay Peninsula-Lingga 


Archipelago and Borneo. This pattern suggests 
that populations in each geographical area might 
be the consequence of post-glacial expansion from 
one or two refugia, but that gene flow between 
different glacial refugia was fairly restricted. The 
present study gives an additional support for 
inference of non-forested land bridge at LGM 
between Sumatra-Malay Peninsula and Borneo. 
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